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Abstract 

I/"") , We develop a constructive method to derive exactly solvable quan- 

' turn mechanical models of rational (Calogero) and trigonometric (Suther- 

■ land) type. This method starts from a linear algebra problem: finding 

eigenvectors of triangular finite matrices. These eigenvectors are tran- 
scribed into eigenfunctions of a selfadjoint Schrodinger operator. We 
^ ■ prove the feasibility of our method by constructing an " AG% model" of 

^ ■ trigonometric type (the rational case was known before from Wolfes 

\ 1975). Applying a Coxeter group analysis we prove its equivalence 

*^ with the model. In order to better understand features of our 

J> ■ construction we exhibit the F4 rational model with our method. 
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1 Introduction 



The completely integrable models are traditionally characterized by their 
relation with simple Lie algebras A n , B n , C n , D n , G 2 , F 4 , E 6 , E 7 , E 8 . This 
relation is the starting point of the Hamiltonian reduction method exploited 
by Olshanetsky and Perelomov These models possess as limiting cases the 
trigonometric (Sutherland) and rational (Calogero) models that are exactly 
soluble, i.e. their eigenvalues and eigenvectors can be derived by elementary 
methods. 

This exact solvability has been shown to follow from the fact that the 
Schrodinger operators can, after a "gauge transformation", be rewritten as 
a quadratic form of Lie algebra operators. These Lie algebra operators are 
represented as differential operators acting on polynomial spaces. This pro- 
gram was formulated in B and successfully applied first to the A n series in 
0. Then it was carried over to the other sequences B n , C n , D n and G2 and 
even to corresponding supersymmetric models ||]. 

Our aim was to turn the arguments around and to develop an algorithm 
which may allow us to construct new exactly soluble models. First investi- 
gations were presented in ||. The program contains two major and separate 
issues, to render a second order differential operator curvature free and to 
find a first order differential operator satisfying an integrability constraint. In 
this paper we present our algorithm in the following version. We start from 
a standard flat Laplacian and introduce Coxeter (or Weyl) group invariants 
as new coordinates. If the Coxeter group contains a symmetric group as sub- 
group, these invariants are built from elementary symmetric polynomials. 
The second order differential operators obtained this way are curvature free 
by construction, and act on polynomial spaces of these Coxeter invariants 
that form a flag. This flag is defined by means of a characteristic vector 
(j?- vector) . 

Then we solve the integrability constraints by constructing "prepoten- 
tials" with a fixed algorithm. These prepotentials define the gauge transfor- 
mation alluded to above which renders the differential operator the form of 
a standard Schrodinger operator of N particles in 1-dimensional space with 
a potential. Each prepotential contributes an additive term to this poten- 
tial with a free (real) coupling constant. Finally the prepotentials define the 
ground state wave function of the Schrodinger operator which originates from 
the trivial polynomial in the flag and thus contains no further information. 
Except a possible oscillator prepotential in the translation invariant cases, 
the prepotentials are in one-to-one relation with the orbits of the Coxeter 
group. 

We show that all known exactly soluble models can be obtained this way 
(at present we have to make an exemption with respect to E 6 , E 7 , E 8 , but 
this will soon be overcome). Applying the method of constructing the Coxeter 
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invariants of A 2 to A3, we obtain an " AG3 model". Its Coxeter diagram 
is that of the affine Coxeter group B 3 , which possesses the same invariants 
as the Coxeter group -B3 . This leads to an explicit proof of the equivalence of 
the AG% model with the B3 model. Thus a translation invariant four-particle 
model after separation of the cm. motion is shown to be equivalent with a 
translation non-invariant three-particle model. In this paper we also discuss 
F 4 from the view point of our algorithm. The Schrodinger operator obtained 
(only the rational case) deviates slightly from the one given in [|IJ (probably 
due to a simple printing error in [p]]). 

Thus our method shifts the centre of interest from the simple Lie alge- 
bras and their homogeneous spaces to the corresponding Weyl groups and by 
generalization to the Coxeter groups. On the other hand, the differential op- 
erators acting on polynomial spaces of Coxeter invariants define Lie algebras 
of their own, but at present these algebras are only of marginal interest. 



2 The constructive program 

We are interested here in the bound state spectrum of Schrodinger operators. 
The whole analysis is therefore performed in real spaces. Consider a flag of 
polynomial spaces Vn(p), N 6 Z>, p e IN n 

V N (p) = span{4 1 z£ 2 ...2£"|npi + r 2 p 2 + ■■■ + r n p n < N} (2.1) 

(Pi e IN) 



We consider differential operators of first order 

dz a 



D?L = (2-2) 



(5 a multi-exponent) 
and of second order 



n(2) _ JSj ° 2 ( 2 Q\ 
D l^- Z dZadzb ^ 



that leave each space Vjy(p) invariant. If 

p= (1, !,...,!) (2.4) 



then the operators ( 2.2|) generate the full linear (inhomogeneous) group of 
1R„ and the operators of second order ( |2.3|) can be obtained as products from 
the first order operators, i.e. in 



a = e (c) , ei c) = 5 c b or a = (2.5) 

and in fl2.3|) 

a = e {c) + e {d) or a = e (c) or 5 = (2.6) 
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Now we consider a candidate for a future Schrodinger operator 

0,c 



D = - E f^a^l^fsw 



',a,b 



The eigenvectors and values of D in Vjv can be calculated easily by finite 
linear algebra methods. Let 

U N = V N /V N -i (2.8) 

and the diagonal part of D on Un be defined as D N 

D N U N = DU N n U N (2.9) 

If the eigenvalues of are all different, the number of eigenvectors equals 
dim?77v But if some eigenvalues coincide (this is true in the generic case!) the 
number of eigenvectors is smaller. Then the Hilbert space on which the final 
selfadjoint Schrodinger operator is acting is not an L 2 -space. The missing 
eigenfunctions can be described. For more details see |J. 

If we want completely integrable models we must make sure that a com- 
plete set of involutive differential operators exists. For this task Lie algebraic 
methods may be very helpful. 

Given a differential operator (|2.7|) one can characterize the vector p in 



2.1 ) by inequalities 



9[S;a,b] ^ =>> pa-p a -Pb<0 (2.10) 

V;c]^° pP~Pc<0 (2.11) 



There should be enough equality signs in (|2.10|) , (|2.11 ) for a chosen p so that 
Dn 7^ 0. It turns out that there exists a minimal p- vector p mm so that the 
^v(Pmin) s P aces are maximal: For each N, p there is N' so that 

V N (p) C MPmin) (2.12) 

It is convenient to work only with this minimal p- vector. 

The first step in transforming D into a Schrodinger operator is to write 
it symmetrically 

= ^-9ab(z)^ + E r.(*)#- (2-13) 

where 

ffi^E^]^ (2-14) 
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We write g~l because this is the inverse of a Riemann tensor. The Riemann 
tensor is assumed to be curvature free. The task to make it so will not 
arise in this work. But we mention that we developed a minimal algorithm 
to solve this issue. 

Following the notations of || we "gauge" the polynomial eigenf unctions 
ip of D by 

if>(z) = e-xM<p(z) (2.15) 

so that 

(g = (det 9 -')-i). 

This is possible if and only if 

E9*(*)-£-[*X-toy/g\=r a (z) (2.17) 

b dZ b 

which implies integrability constraints on the functions {r a (z)}. If they are 
fulfilled we obtain a "prepotential" 

p = lnP (2.18) 

so that 

p = 2 X -ln v ^ (2.19) 
In most cases studied, we found solutions for p as follows. Let 



det g- 1 (z) = 11 P, t (z) (2.20) 

i=l 

where {Pi(z)} are different real polynomials. Then 

p(z)=j2li^Pi^) (2-21) 



i=l 



with free parameters ji solves the requirement that {r a (z)} ( [2. 17[ ) belong to 
differential operators leaving each Vn invariant. In particular 



-(0 



1 ^ ,dR 



rx'(z) 



am (2 - 22) 

are polynomials. Inserting fl2.2U| ), ( |2.21| ) in ( 2.19 ) we obtain finally 

X=\jl{li-\)^Pi (2-23) 
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We will later see that in the case of the models of Calogero type a term 

7ohiP (2.24) 

can be added to p, where 

P {z) = e 21 (2.25) 

is not contained in detg -1 as a factor. This prepotential gives rise to the 
oscillator potential. 

Finally we mention that e~ x is the ground state wave function of the 
Schrodinger operator, as follows from (|2.15|) . 

The expression ||, (6.17) for the potential W(z) contains a term linear 
in x 

Each divergence 

H4-rf{z) = C^ (2.27) 

ought to be a constant. From now on we shall dismiss all constant terms in 
W(z). 

We can then write the potential as 



W(z)=J2^ R n(z) (2-28) 

(2.29) 



Kij ^ 9ab dz a dz b 



lii = \(lili-\) (hJy^O). (2.30) 
In the cases of this article 

Rij = const if i ^ j (2.31) 



If we then set 



we obtain 



H = -Vi + \ (i^O) (2.32) 



with 



W(z) = ^-y ii R H (z) (2.33) 
i=i 

7« = ~ 1) (2-34) 

As stated in the Introduction the variables {z{\ appearing in this section 
are identified with Coxeter invariants formed from root space coordinates 
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{x n } or {y n }. These invariants are either polynomial or trigonometric. Fi- 
nally we return from the invariant coordinates {z,j} to the root space coordi- 
nates {x n } in the Schrodinger operator (|2.16|) . Each contribution 

Rii = % (2.35) 

admits a partial fraction decomposition due to the factorization of the prepo- 
tentials Pj (Section 5). The label i — 1 is always reserved to a "Vandermonde 
prepotential" , i.e. 

Pi ~ Y\.( x i ~~ x i) 2 or XT( s ^ n ( a:i ~~ x i)) 2 (2.36) 

i<j i<j 

or alike. 



3 Translation invariant models 

3.1 Relative coordinates 

The Laplacian for an Euclidean space Hjy 

N Q2 

A = T.£t (3-D 

i=l ux i 

is translation invariant. We introduce relative coordinates by 

Vi = ^~Jj X ( 3 - 2 ) 

JV 

X = J>< (3.3) 

i=l 



They separate the Laplacian such that 



d 2 N d 2 1 / N d \ 2 
A = ^ + E^-4fe^ (3-4) 



We use all {yi}f =l as coordinates on the plane 



N 

E^ = ° (3-5) 
t=i 



in order to maintain permutation symmetry. 
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3.2 Elementary symmetric polynomials 

Elementary symmetric polynomials of N variables {qi}f =1 are defined by a 
generating function 

N N 

J2Pn(q)t n = U(l + qit) (3.6) 

n=0 1=1 

They are invariant under the symmetric group Sjy. For each g e Sn we have 
a sector (simplex) E g C Hat 

^9 = {fti < fe < • • • < Qi N ; in = g{n)} (3.7) 

so that 

IRtv= [j E g (3.8) 
Inside we can use the {p n }n=i as coordinates since 

M m = ^ (3.9) 

r JV i 

detM = (-l) [ - ] V( qi ,q 2 ,...q N ) (3.10) 
where 1/ is the Vandermonde determinant. 

3.3 The An-i series 

The root system of A N _ X and the corresponding Weyl group possess elemen- 
tary symmetric polynomials as invariants. We express the Laplacian in each 
sector Eg ( p.7|) intersected with the plane ( |3.5|) in terms of these polynomials 



Tn(yi,-,VN) =P»(9)|«=WiiJli (3.H) 

The dynamics will be bounded to such sectors by corresponding potential 
walls automatically. 

Then (see [3]) it results 

v — -1/vAV 

N q2 N q 

= E^ + E^ (3-12) 

n,m=2 u! nVl m n= 2 u 'n 



with 



and 



9nM = ^(m -l)(N-n+ l)r n r m - T n ^ m ^(r) (3.13) 

T nm {r) = ^2(21 + n- m)T n+l T m _i (3.14) 
i>i 
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Here it is understood that 



T = 1 

n = o 

T n = Oforn<0, n>N (3.15) 

In this case det g~ x is indecomposable as a polynomial, so we set 

P = e UT2 (3.16) 
P l = det^T 1 = C N V( yi ,...,y N ) 2 (3.17) 

The resulting vectors {r a }^ are 

r<® = (-2r 2 ,-3r 3 ,...,-Nr N ) (3.18) 

: explicit formulas known only for N < 4 (3.19) 

and the potential is 

1 1 N 

-W(x) = -u^^ + 9 E (Xi-*i)- 2 (3-20) 

Z Z i=l l<i<j<N 

The corresponding Sutherland models are obtained as follows. We use as 
coordinates a system {o"„}^ =2 defined by (these differ from those in [3]) 

N 

o"o = ri cos ^ ( 3 - 21 ) 



i=i 



and 

The identity 



CTn = 00 • Pn(5)|«=tan W (3.22) 



1 = exp ^ E % 

AT 

= n( cos %' + isin %) 

AT 

= J2^n(y) (3.23) 

n=0 

allows us to eliminate <To and o\ in terms of the remaining {(J n }n=2 so that 
polynomials go into polynomials. 
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The Laplacian is expressed correspondingly as 

N Q2 x / N q \ 2 

at w 5 

= ^ 9nm d^~do~ + ^ ^"do 7 

n,m=2 UU n UU m n=2 UU n 



9nm( a ) — —Tn+l,m+l(&) — T n+ l !m -l(a) 
— ^n-l,m+l(°") — ^~n-l,m-l(°") 

+ ^[(m + l)cr m+ i + (m - l)<r m _i] 

x [(at - n - l)a n+1 + (N - n + lK_ x ] (3.25) 

with T nm as in ( J3.14Q . 

Once again det g" 1 is indecomposable, so we set 

P 1 = detg- 1 = C' N V(y 1} ...,y N ) 2 (3.26) 

where 

V(y h ...,y N ) = Y[sin(yi - yj) (3.27) 

i<j 

has the symmetry of the Vandermonde determinant (translations and per- 
mutations). The vector is known only up to N = 4. Finally we obtain 
as potential 

\w{x)=g J2 sin(x, - x,)- 2 (3.28) 

l<i<j<N 

In each case Apj-i the minimal p-vector is (1, 1, 1) G IN^ -1 . 



3.4 The G 2 and AG 3 models 

The models G 2 and AG% belong also to the domain of translation invariant 

models [§J. For G 2 we start from A 2 and extend its Weyl group by a Z 2 
group 

Vi -> 

As invariant variables we use j| 

A 2 = r 2 (3.29) 

A 3 = r 3 2 (3.30) 
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In these variables 



We find 



so that 



kdy? Ah dm) 

3 Q2 3 q 

z X] Sab a, q> + XI 



-2A 2 , -6A 3 

-e 

4 



— 6A 3 , +IA2A3 



det( ? - 1 = --A 3 (4A3 + 27A 3 ) 
Thus as ansatz for the prepotentials we use 

P = e wX2 

Pi = 4A^ + 27A 3 

P2 = A3 

The r- vectors (justifying this ansatz) are 

r (0) = (-2A 2 ,-6A 3 ) 
r M = (-6,0) 

r {2) = (-6,+^) 

The minimal p-vector is 

P=(l,2) 

The potential is 

i=i 

l<i<i<3 i<j,W(i,i) 

with 

01 = ^(fl - 1) 

5-2 = 3z/ 2 (z/ 2 - 1) 

If 

z/ 2 = or u 2 = 1 

we return to the A 2 model. 
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In the Sutherland case we use as variables 



\ii = o 2 (3.44) 
/is = °l (3.45) 



leading to the inverse Riemann tensor 



„-i _ ( -2^2 - 2$ + §M3, -Ma(6 + f M2) \ ^ 

5 -V-^3(6 + ^ a ), 1^3-8^ J ^ 

Now det g -1 is decomposable with 

detg- 1 = -^P 1 (^) (3.47) 

and 

Pi(/i) = 4/x§ + // 3 (8/x| + 36/i 2 + 27) + 4/4(1 + /i 2 ) (3.48) 

P 2 ( M )=/i 3 (3.49) 



The r-vectors are 



r (1) = (-6-8/^,-16^) (3.50) 
16 8 
Y^ 2 ' 3' 



16 8 

r (2) = (-6 - —// a , -axI - I6//3) (3-51) 



The resulting potential is 
1 



2 



W(x) = gi sin(xj — .r 



2 

i<«<i<3 

+ XI sin ^(xj + - 2x fe ) -2 (3.52) 

In the case of the A 2 models the spaces Vjv decompose into even and 
odd subspaces in r 3 (or a 3 ) which are left invariant separately under action 
of the Laplacian. In the case of the odd spaces we can factor ^(a^) and 
leave an even space as well. In each case we obtain a polynomial space in 
the variables A 2 , A3 = r| (/x 2 ,/i3 = of). Thus starting from such polynomial 
space and multiplying with T3 2 (o^ 2 ) we obtain the A 2 model if z/ 2 = or 
z/ 2 = 1 but a new potential in all other cases. 

It is plausible that a similar procedure works for A% but not for A/v-i, N > 
5. In the latter models we have two or more odd variables r 3 , T5, ...(03, 05, ...) 
and there is no factorization of the odd invariant subspaces. Let us sketch 
the A 3 model whose extension leads to the AG% model ||. 

In this case the variables are chosen as in (3.29), (3.30), (3.44), (3.45) 

A 2 = r 2 , A 3 = r 3 2 , A 4 = r 4 (3.53) 



12 



The inverse Riemann tensor is 
/ -2A 2 , 



9 



-6A 3 , 



-4A 4 

-6A3, 4A 3 (A 2 ! -4A 4 ), A 2 A 3 
V -4A 4) +A2A3, -2A 2 A 4 + |A a 



The determinant is decomposable as 

det^- 1 = A 3 Pi(A) 
and the ansatz for the prepotentials is 



Po(A) = 
A(A) = 

P 2 (A) = 

The r-vectors come out as 

r (o) 
r (i) 

~(2) 



27\ 2 3 - 256\l + 128A^A^ 
-16A^A 4 + 4A^A 3 - 144A 2 A 3 A 4 
A 3 



;-2A 2 ,-6A 3 ,-4A 4 ) 

;-12,0,-2A 2 ) 

;-6,4(A2-4A 4 ),A 2 ) 



The potential for this Calogero type model is 



1 



W(x) 



-w 



E 



x- 



1=1 

+91 E < 

l<i<j<4 



+ 92 E ( x i + x i ~ x fc ~ 



3 terms 



with 



9i = v\iy\ - 1), 92 = 2v 2 {v 2 - 1) 

It was discovered first by Wolfes, [[F]. 

The Sutherland model is obtained in the same fashion. With 

the inverse Riemann tensor is 

g 2 2 = — 2yU2 - 2/^2 - 8/X4 + 2;U 3 + 8;U 2/ u 4 + 8/^ 

fe 1 = -6^3 - 4/i 2 /i 3 

^ = -4/i 4 - 6/i 2 /i 4 + /i 3 + 4/i4 

9 3 3 = 4/i 3 [-4/i 4 + /i2-4/i 2 /i 4 + 4^-2/i 3 ] 

fi^ 1 = V2H3 ~ 6/x 3 /i 4 
-1 3 

9U = - 2 ^2/i4 + ^3 
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Its determinant decomposes 

det^T 1 = -frPM (3.71) 

Pi(fj) = 256/4 + 32 further terms (3.72) 
(equ. (A. 2) from [g]) 

P 2 (/i) = /is (3.73) 

and the r-vectors are 

r W = (_l6/i 2 - 12, -24/x 3 , -12/^4 - 2/ia) (3.74) 

r*- 2 - 1 = (— 4/x 2 — 8, 16/X4 — I6/14/X2 + 4/i 2 — 8/13 — I6//4, — 6/X4 + /i 2 ) (3.75) 
The factorization of 03 which is necessary in this case is 

0-3 = - [J sm(^ + yj) (3.76) 

l<i<j<3 



implying 

(hl -4 2 (sin(^ + 2/j ))- 2 (3.77) 



P 

2 i<*<i<3 
This gives the potential 

-W(x) = g x Yl (sm(xi- Xj)y 2 

l<i<j<4 

+ 4#2 (^-(xi + Xj - x k - Xi))~ 2 (3.78) 

3 cases 

The discussion of this model is resumed in Section 5. 



4 Translation non-invariant models 

4.1 The BC N and models 

As we shall see there is only one series with two (Calogero) and three (Suther- 
land) independent coupling constants. For any such model we use as Carte- 
sian coordinates {xi}f =l and require permutation symmetry Sn and reflection 
symmetry {~%-2) N x i — ► ~ x % for each i separately. Then the natural coordi- 
nates invariant under these group actions are || 

An (a?) =Pn(g)|«=»J,aUt ( 4 - X ) 

There is a bilinear relation with the {p n {%)}n=i 

2n 

K(X) = Y,(- 1 T' k P2n-k(x) Pk (x) (4.2) 
k=0 
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The inverse Riemann tensor for the full Laplacian (|3.1|) is then 

g-^(X) = 4M„ m (A) (4.3) 

where we introduce the shorthand 

M nm (\) = X](2/ + n - m + l)A n+i A m -i-z (4.4) 
z>o 

Its determinant factorizes 

det g- 1 = (-lp]A N X N P 1 (X) (4.5) 

where 

P 1 (X) = N N \^- 1 + ... (4.6) 
= D N V(xl,xl, ...x 2 N ) 2 

and 



P 2 (A) = A^ (4.7) 

Both functions Pi,P% factorize in a trivial way. In the general case there is 
no explicit expression for A 1 ' but 

= 4(iV-a + l)A a _i (4.8) 

If follows 

^ 22 = 4^ = 45X 2 (4-9) 
Av i=l 

The resulting potential is, including an oscillator potential 
1 1 N 

-W(x) = -u 2 Y.^ + 9i E [(Xi-x^ + ixi + xj)- 2 ] 

Z Z i=l l<i<j<N 

N 

+^E^r 2 (4.io) 

i=l 

<?! = z/i(z/i-l) (4.11) 
3 2 = ^2(^2-1) (4.12) 

In the Sutherland case we use coordinates 

N 

fi = Y[ cos 2 Xi (4.13) 
i=i 

/i«(x) = /^(^Pn^lg^tan^.aiH (4.14) 

n G {1,2, ...N} 
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From the identity 



N 

1 = I^[(cos 2 Xi + sin 2 Xi) 



i=i 

N 



n=0 

we learn how to eliminate /x in facour of {/in}^ so that a polynomial of 
{f^n}n=o remains a polynomial. 

In this case the inverse Riemannian is 

9nrn = 4{M n+1>m+1 (fj.) + M n>m (fl) 
,m+l 

(/i)-M n+1 , m (/i)} (4.16) 

and the determinant decomposes as 

det g- 1 = A N (-lp]fi of i N P 1 (fi) (4.17) 
Now the factorization of P\ (/x) is 

Pi(fi) — D' N Yl (cos 2 Xisin 2 Xj — sin 2 Xicos 2 xj) 2 (4.18) 

l<i<j<N 

and we choose 

P 2 (/i) = Hn (4.19) 

P 3 (/i) = Mo (4.20) 

Again we have no general explicit expression for but 

r< 2 > = 4[(7V-a + lK_ 1 -(7V-aK] (4.21) 

= 4[(a + lK +1 -a// a ] (4.22) 

so that 

R 22 = ^zi = 4 Ecot 2 a; J (4.23) 

i=1 

U ^ 

#33 = — = 4Etan 2 a; l (4.24) 

Thus we end up with a potential 

^W(a:) = <7i E [( sin (^ ~ x j))~ 2 + ( sin (^ + ^))~ 2 ] 
i<i<j<N 



JV 



+22 E 1 



sin 



i-2 



1=1 

N 



+g 3 J2( C0S ^y 2 ( 4 - 25 ) 
1=1 
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where gi >2 are as in ( |4.11| ),( f4~T2|) and 



93 = ~ 1) (4.26) 

An alternative form of the potential is obtained from 

92 93 92 ~ 93 4:93 97 x 

sin x cos z x sin x sin 2x 

If we set g>2 = #3 or g 3 = we obtain different samples of the BCn or Dn 
series. We mention finally that the minimal p-vector is in all cases 

p = (1, 1, ...1) G JN N (4.28) 

4.2 The F 4 model 

The F4 model belongs also to the translation noninvariant class. The Weyl 
group of F4 possesses four basic polynomial invariants 

h(x), I 3 (x), h(x),I 6 (x) (4.29) 

(I n of degree 2n) which can be expressed as polynomials in the {X n } n=1 as 
follows 

(4.30) 

^AiA 2 (4.31) 

— Ai A 3 -|- — Ao 
4 3 12 2 



h 


= Ai 


h 


= A 3 - 


h 


= A 4 - 


h 


= A4A2 



1 \3 , 1 \2\2 1 



(4.32) 
(4.33) 



36 z 24 z 1 64 1 
In these coordinates the inverse Riemannian can be given as 

g£ = Aml m (4.34) 

933 1 = fA/i-^si? (4-35) 

1 3 3 

fo 1 = 8I 6 -3li--Ijf--hll (4.36) 

g£ = 16/ 4 2 + + 14 Wi + \lll 2 x - \hl\ - (4.37) 

9ll = -4/ 4 /3-2V 1 + ^ 4 / 1 3 + ^/ 1 + i 3 / 1 4 (4.38) 

4 4 16 

9* = Mlh + 2hhll- l -hll (4.39) 
g£ = 30 Wi + yWi - ^Vf + 12/ 4 2 / 3 + 6/4/3^1 
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The determinant decomposes into two factors 

d6t ^ 1 = 3^72 Pl(/)jP2(/) (441) 

where Pi (I) is connected with the Vandermonde determinant squared as 
usual 

Pi (J) = -4096/ 4 3 + 432/ 4 + 3072/ 2 - 2304/ 6 / 4 / 2 
-hlQhhll + 864/ 4 / 3 2 / 2 + 2I6/4V1 5 

+432/ 2 / 4 + 27/lit - 2304/ 6 / 3 2 + 216/ 3 3 / 3 (4.42) 

or in factorized form 

Pi(/) = -16 LI (^ 2 -^ 2 ) 2 (4-43) 

l<«<i<4 

and P 2 (J) 

P 2 (I) = 36864J 2 - 18432 hUl\ - 4608 V3 A 3 + 32I 6 lf 
-49152/ 3 - 36864J 4 2 / 3 / 1 + 1536/ 2 / 4 
+768/4/3/1 - 12/ 4/i 8 - 9216/4/3/1 
-768/|/i 3 + 96/ 3 2 /i 6 - 3/3/1 9 (4.44) 

which factorizes as 

P 2 (/) = -12A 4 (64A 4 -16A 2 : + 8A 2 A 2 -A 4 ) 2 

4 

= -12x 2 x 2 x 2 x 2 J] (^i-El- 1 )" 1 ^) 2 (4-45) 

The r-vectors are 

r« = (48, -211 0, 36/ 4 /i + I2/3A 2 - ^II) (4.46) 

16 

r (2) = ( 48) _ 4 /2 ) _i2/ 3 , 24/4/1 + 6/ 2 / 3 - (4.47) 

o 

The potential resulting is 

-w{x) = -u 2 x l + 9i E [(^-^r 2 + (^ + ^r 2 ] 

1<«<4 l<i<i<4 

+£ 2 { £ +Ex" 2 } (4.48) 



e{+i,o} 



i=2 / i=l 



where gi j2 are as in ( |4.11| ),( f4.12|) . The minimal p- vector is 

P= (1,2,3,5) (4.49) 



18 



5 Coxeter groups, orbits and prepotentials 

The prepotentials used in the empirical constructions of sections 3 and 4 ne- 
cessitate a mathematical interpretation. Let W be a Coxeter group generated 
by the reflections 

{s a } (5.1) 
where a are roots running over a set 

$ = {«}f (5.2) 

The roots span an Euclidian space V. In this space the reflections {s a } act 
by 

x G V : s a x = x — 2 7^ — -a (5.3) 
(a, a) 

If the Coxeter group W is " crystallographic" , it is a Weyl group (for more 
details see ||). 

We denote a set of basic polynomial invariants of W by 

j>i(x), . . . ,z n (x)}, n = dimV (5.4) 

Invariance means 

Zi(w~ l x) = z,i(x) 

= wzi(x) (5.5) 

for all w G W . The Jacobian for the transition {xj} — ► {zj} 

^det{||} (5.6) 

can be factorized as follows ([[J, Proposition 3.13). 

Each reflection s a leaves a hyperplane H a in 1/ pointwise fixed, let H a be 
given by a linear function Z Q 

l a (x) = (5.7) 

Then due to the proposition 

J = C I] (5-8) 

with $ + the set of positive roots. The proof of this proposition is rather 
elementary. 

For any inverse Riemann tensor {g~ 1 } of Sections 3 and 4 we obtain this 

way 

det^ = C 2 I] W (5-9) 
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If $ decomposes into orbits under W 

$ = U$i ( 5 - 10 ) 

i 

then 

Pi = n ( 5 - n ) 

is an invariant polynomial under action of W and therefore a polynomial in 
the basic invariants 

P i = P i (z l ,...,z n ) (5.12) 

These polynomials are the prepotentials constructed in Sections 3 and 4. The 
factorization of these prepotentials as quoted at the end of Section 2 (eqns. 
(|2.35|) , (|2.36|) ) and used throughout in Sections 3 and 4 is based on ( |5.11| ). 



We emphasize that our empirical results of Sections 3 and 4 indicate the 
validity of further mathematical propositions which could not be traced in 
the literature: 

1. an analogous factorization theorem for the trigonometric invariants; 

2. the polynomial properties ("integrability") of the functions r^(z) ( |2.22| ). 

Now we return to the AG% model of Section 3. We identify the roots 
involved in a model using (|5.7|) , (|5.9|) 

l a (x) = (a v ,x) 
In 

(a v = -, the "dual" of a) (5.13) 

(a, a) 

and the Sutherland version whose potential is 

\w{x)= £ 9i £ [sin/ a (x)r 2 (5.14) 
orbits i a ^t 

Thus the simple roots of A 3 



a i 


= ei 


- e 2 






= e 2 


- e 3 


(5.15) 


a 3 


= e 3 


- e 4 





are completed by a fourth root in AG% 

a A = e 3 + e 4 - ei - e 2 (5.16) 
The corresponding Coxeter-diagram is shown in Fig. 1. It belongs to the 



20 



1 




Figure 1: Coxeter diagram of B 3 

affine Coxeter group B 3 (||, Figure 1 in Section 2.4). 

The coordinates of the B 3 root space with respect to the standard basis 
{fi}i=i are denoted {£i}f =1 , those of AG 3 with respect to the standard basis 
{ e i}t=i by {%i}f = i as before. The simple roots of B 3 are 

Pi = fi-h, Pi = h~h, fo = h (5.17) 
and ^3 is obtained by adjoining 

Pa = ~fi ~ H (5.18) 

It follows that 

sa 16 I = I -6 I (5.19) 




£3 

leaves the Coxeter invariants of B 3 

Ai(0 = E *? (5-20) 



l<i<3 

a 2 (0 = E e'e? (5-2i) 

i<«<i<3 

As(0 = (5-22) 

invariant, too. This suggests the equivalence of the AG 3 and the B 3 models. 
An explicit identification of the simple roots 

fi = ^(ei - e 2 - e 3 + e 4 ) (5.23) 

/ 2 = ^(-ei + e 2 -e 3 + e 4 ) (5.24) 

/s = \{~ei -e 2 + e 3 + e 4 ) (5.25) 



gives (i, j e {1,2,3}) 



£i - aJj = & - (5.26) 

x 4 - % = E & ( 5 - 27 ) 
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A 2 G 2 G 2 

adding a reduction 
reflection 



Figure 2: Extending the Coxeter diagram of A 2 to G 2 and reduction to G 2 
It follows 

l<i<j<A 3 cases ^ 

= 9l £ {[sinte-^r 2 + [smte + ^r 2 } + ^ 2 E[siner 2 

1<*<7<3 i=l 



(5.28) 



Moreover the rational invariants ( ft.64|) can be identified with the invariants 
(|5T20|)- (|5T2l) 



li 2 (x) = ~Xi(0 (5.29) 
l*(x) = (5.30) 
N( x ) = -\\ 2 (0 + ^Xi(0 2 (5-31) 

This establishes the equivalence between the two models. 

Our method involves a reduction of the affine Coxeter group B% to the 
Coxeter group B 3 having the same invariants. It may therefore be of interest 
that the construction performed in M is analogous (see Fig. 2). 
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